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ended,  the  "cross  section"  is  shown  to  fall  off  as  the  recipro¬ 
cal  of  the  time  squared. 

Finally,  an  estimate  is  made  of  the  effect  of  a  non-uniform 
electron  density  and  it  is  found  that,  provided  the  rate  of 
decay  is  not  too  great,  the  cross-polarizcd  backscatter  and 
the  directly-polarized  backscatter  arise  from  the  same  volume 
of  the  scatterer. 
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TIME  BEHAVIOR  OF  MULTIPLE  SCATTERING 


SECTION  1 

INTRODUCTION 

Ruffine  and  de  Wolf*  have  calculated  the  steady-state,  incoherent,  cross-polar¬ 
ized  scattering  of  electromagnetic  waves  by  a  turbulent,  tenuous  plasma  using 
the  second  Born  approximation.  The  principal  result  of  this  calculation,  under 
the  conditions  that  these  authors  assume,  is  that  the  main  contribution  to  the 
cross  section  for  this  process  comes  from  terms  corresponding  to  the  first  and 
second  scatterings  occurring  in  different  correlation  cells.  Thus  the  effect  of 
statistical  averaging  is  solely  that  of  determining  the  magnitude  of  the  scattered 
amplitude  from  each  cell. 

The  usefulness  of  cross-polarized  data  in  the  analysis  of  radar  returns  from 
turbulent  media  is  limited  unless  it  can  be  shown  that  most  of  the  second  scat¬ 
terings  take  place  within  some  maximum  distance  of  the  first  scatterings.  If 
they  do  not,  then  the  returned  signal  will  be  very  much  stretched  out  in  time. 

The  signal  will  have  "sampled"  a  large  portion  of  the  plasma  and  one  could  not 
then  assign  a  localized  cross  section  for  cross-polarized  scattering. 

In  this  report  wc  consider  an  idealized  form  for  the  wake,  namely,  a  circular 
cylinder  of  diameter,  d.  As  is  readily  seen,  even  directly  scattered  signals,  sup¬ 
posedly  given  accurately  by  the  first  Born  approximation,  will  be  stretched  out  by 
a  "time  stretch"  of  the  order  2  d/sin  0,  where  0  is  the  angle  between  the  direction 
of  propagation  of  the  incident  wave  and  the  axis  of  the  cylinder.  We  neglect  of 

I)  R.S.  Ruffine  and  D.  A.  de  Wolf,  Cross  Polarized  Electromagnetic  Backscat- 
ter - hereinafter  referred  to  as  (1 ) 
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course  the  very  small  sphericality  of  the  incident  wave  at  the  position  of  the 
plasma. 

In  this  report  we  also  use  Hie  second  Born  approximation  but  consider  an  input 
of  the  form  of  a  square-wave-modulated  r-f  signal.  The  carrier  frequency  is 
taken  to  be  higher  than  the  plasma  frequency  of  the  scattering  medium.  The 
time  dependence  of  the  scattered  power  is  investigated  for  times  which  corres¬ 
pond  to  pulse  stretching  many  times  longer  than  the  unavoidable  stretch  men¬ 
tioned  in  the  previous  paragraph.  A  comparison  is  made  between  the  power  re¬ 
turned  at  such  times  and  that  returned  at  earlier  times.  The  principal  result  is 
that  the  "cross  section"  decreases  as  the  inverse  square  of  the  stretch  time,  for 
long  enough  times.  Also  the  cross  section  becomes  increasingly  aspect  sensi¬ 
tive  as  time  increases. 

In  Section  II  we  derive  the  relevant  field  amplitude  and  discuss  some  general 
aspects  of  the  problem.  In  Section  III  wc  formally  introduce  the  statistical  aver¬ 
aging.  The  so-called  isochronal  surfaces  are  discussed  in  Section  IV.  In  Sec¬ 
tion  V  we  compute  the  incoherent  cross  section  for  large  stretch  times  and  com¬ 
pare  the  results  with  those  of  (1).  Finally,  in  Section  VI  we  estimate  the  proportion 
of  the  scattering  which  arises  from  the  resolution  cells  to  the  scattering  from  one 
resolution  cell  when  the  electron  density  is  not  constant. 

In  Figure  1  there  is  shown  a  schematic  of  the  scattering  geometry  and  the  rela¬ 
tive  sizes  of  the  various  distances  involved  in  the  problem. 


Figure  1.  Schematic  of  the  Scattering 
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SECTION  2 


THE  WAVE  EQUATION  AND  ITS  SOLUTION 


We  use  the  microscopic  form  of  Maxwell's  equations  in  Gaussian  units 


Vx  E  - - H 

c 


_  — -  1  9E 

V  x  H  =  —  3  +  —  r— 
C  C  ot 


(1) 


For  a  tenuous  plasma  we  use  the  constitutive  equation 


ii 

9 1  dt 


2  - 
co  E 

p 


(2) 


2  2 

where  C0p  =  4  77  n  e  /m  is  the  plasma  frequency  (n  is  the  electron  number  den¬ 
sity.  )  Following  the  usual  procedure  of  eliminating  the  magnetic  field  one  gets: 

,2 


1  9  -=*■ 

—  —  E+VxVxE 

c2  dt2 


CO 


p 


(3) 


We  wish  to  solve  this  equation  for  the  electric  field  by  converting  it  into  an  in¬ 
tegral  equation  and  by  then  using  the  first  few  terms  of  the  corresponding 
Neumann- Liouville  series  (successive  Born  approximations. )  In  order  to  do  this 
it  is  necessary  to  invert  the  vector  operators  on  the  left-hand  side.  We  shall  do 
this  symbolically.  For  didactic  reasons  we  shall  first  treat  die  steady-state 
equation  considered  in  (1).  This  equation  can  be  written 

-VxVxE  +  k2E=f  (4) 

If  E  were  purely  transverse,  then  the  solution  would  be 
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E  = 


V2  +  k2 


f  j  j  /  j  /  3  /  _2  ,  2  / 

f  =  J  G(r,  r  f(r  )  d  r  (V  +  k  )  G(r,  r  ) 


2  ,  2 


=  (V*  +k")  G  (r,  r')  =  6(3)  (r  -  r') 


(3) 


As  both  E  and  f  possess  longitudinal  parts  we  first  decompose  equation  (4)  into 
such  parts  and  get 

-*  1  ->  fL 

ET  =  — - -  et,  =— r  (b) 


V2  +  k2 


This  follows  since  x  v  x  E  is  purely  transverse  and  also  since 


V  (V  •  Et)  =  0. 


Now 


f  =  v  _L  v'  • 


f 


M—L. 

J\4TT  V- 


r 


V'  •  fir')  d3  r' 


(?) 


Combining  the  results  of  equations  (6)  and  (7)  one  gets 

f 


E  = 


1  t  L 

_  fm  +  - 


V2  +  k2  “  k*  +  k 

This  last  equation  can  be  written 

1 


i — |"f-v-  v'-fl^ 

+  k2L  V2  j  k2  A2 


9  I  /  1  1  /  r 

-  V  v  -  V  .  f  =  — -  V  — - -  v  •  f 

y2+k2  k2(k2+v2)  £  k  k  +  A 


1  f  = 


The  last  term  stands  for 


-  vf  G  (r ,  r /)  V'  *  f(r  ')  d3  r  '  -  V  J(  ’  •  G)  fd3  r' 


(8) 


(9) 


—  V  V'  • 

i  2 


_2  2 
V  +  k 


(10) 


The  first  equality  follows  upon  integration  by  parts  and  the  parenthesis  indicates 
that  the  differentiation  is  to  be  performed  only  on  die  Green's  function  G.  In¬ 
serting  the  results  of  (10)  into  (9)  yields  (1  -  (5)). 


Returning  to  the  ease  at  hand,  we  have  an  equation  of  the  form: 

-VxVxE-D2E  =  f, 


P  1 

D2  =  -  — 

e2  St2 


(11) 


c ov 
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Proceeding  as  before  we  have 


E-p  = 


- - - £p;  E,  =  -  ~  ■  t 

-  Dz  D 


(12) 


These  equations  are  quite  similar  to  those  in  (6);  the  operator  D  has  replaced 
2  2 

-k  .  Sinec  D  commutes  with  V“,  we  ean  immediately  write 


E  = 


A  - 

\  D“  /  V2  -  D" 


or  equivalently 


E  (r,  t) 


(-^)  / 


G  (x,  x  ')  f  (x  ')  d^  x  ' 


where 


4  3 

x  =  (r,  et),  d  x  =  d  redt  and  G  is  given  by 
G(x,x')  =  47r|i__/|6^  -r'l  -e(t-t')J,  t>t' 


=  0 


,  t  <  t 


(13) 


(14) 
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The  operator  D-2  js  an  integral  operator  inverse  to  D2. 


Denoting  the  operator  1  -  — ^  by  P  we  write  down  the  first  three  terms  of  the 
familiar  Born  approximation: 


E(x) 


E0  (x)  +  4tt  r0  P  (x)  •  /g  (x  -  y)  n  (y)  E0  (y)  d4  y 
+  (4  77rQ)2  P  (x)  J  G  (x  -  y)  n  (y)  P  (y)  J G0  (y  -  z)  n  (z)  d4zd4y 


ro 


m  c 


The  third  term  in  the  above  expression  gives  the  lowest  order  contribution  to  the 
eross-polarized  field  at  far  distances  from  the  seatterer. 


The  eross-polarized  amplitude,  in  the  seeond  Born  approximation,  at  the  posi¬ 
tion  of  the  receiver  is  then 


ee  1  E 


where 


=  £  =  (4Trr0)2jG(x-y)n 


00 


A  A 

(ee  -V)  (e0 


V) 


G  (y  -  z)  n  (z)  E0  (z)d4 


i4y 

(16) 


and  eQ,  ee  are  the  unit  polarization  veetors  for  the  clireet-  and  eross-polarized 
directions,  respectively.  We  have  assumed  that  the  seatterer  is  far  from  the 
receiver  and  have,  thus,  retained  only  the  transverse  part  of  the  scattered  field. 


In  what  follows,  we  shall  consider  that  the  transmitter  and  receiver  are  loeated 
at  the  same  point  and  we  shall  take  this  position  to  lie  at  the  origin  of  the 
(spherical)  eo-ordinate  system.  The  simplest  next  step  that  suggests  itself  is 
to  consider  an  incident  field  of  the  form 

E  =  eQ  E0  -  eQ  |-6  (t  -  r/e)  (17) 
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and  then  to  find  the  resultant  scattered  field  as  a  function  of  time.  However,  such 
an  input,  if  we  use  the  Born  formalism,  would  result  in  the  plasma  particles 
seeming  to  gain  a  constant  speed  and  undergoing  an  unbounded  displacement  in¬ 
creasing  linearly  with  time.  Now  an  actual  plasma  subjected  to  such  an  electric 
field  would  undergo  (damped)  plasma  oscillations.  These  oscillations  are  the 
resonances  of  the  system;  it  is  well  known  that  resonance  behavior  cannot  be 
described  by  any  finite  number  of  Born  terms.  Instead,  we  must  consider  in¬ 
puts  that  are  higher  in  frequency  than  the  plasma  oscillations.  For  this  reason 
we  take  an  input  of  the  form 


E0(r,t)  - 


e  f 
co 1 

R 

A  f 


sin  co 


(  r\ 

(  r  \  /  r  \ 

It  —  1 

H  lt  +  T--  -  H  It - 1 

\  c/ 

L  \  c  /  \  c  /  _ 

(18) 


S  1  ;  H<u)  = 


u  >  0 
u  <  0 


i.  c.  ,  a  square-wave-modulated  pulse  of  r-f  at  the  angular  frequency,  CO.  At  the 
position  of  Hie  radar  set,  r  =  0,  the  pulse  begins  at  t  =  -  T  and  ends  at  t  =  0. 

We  choose  T  such  that  CO  T  is  a  multiple  of  2  77. 


The  amplitude  at  the  receiver  can  now  be  written  more  explicitly  after  integrat¬ 
ing  by  parts  as 


R 


f fd  l’ 2dCt2  6  [r°  “  C  (t~  "  to) 

f 


n  (r  ) 
2 


(19) 


A  A  „ 

6t  •  V2  6o  '  V2 
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12 


■»1»  -  1  J  rl\ 


d  r ,  c  d  t 
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Here,  R12  =  |r.  -  r2  |  and  R  is  the  (almost  constant)  distance  from  the  receiver 


to  a  point  on  the  cylinder.  Carrying  out  the  time  integrations,  wc  get 
2 


5  “R2(/d3r2n<r2>  V  V2  V  V» 


J  Rl2  d2  \  c  c  c  /  1 


(20) 
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It  is  to  be  noted  that  the  differential  operators  act  only  on  the  function  R 


12' 


Beeause  of  the  nature  of  the  function  I,  —  I  will  be  zero  for 

D2 


R- 


T  +  t  -  —  -  — 
e  e 


<  0, 


a  function  of  time  for 

rl  +  r2  +  R12 

-  T  + - 


<  t  < 


rl  +  r2  +  R12 


e  e 

and  a  constant  for  t  >  (r-^  +  ro  +  R^/c.  Recall  that 

a  Jo  32  po 


Wp  E0  "  g l 


at2 


(21) 


(P0  -  eleetrie  polarization)  so  that  measures  the  polarization  of  the  medium. 
Ill  the  Born  formalism,  then,  there  exists  a  remanent  polarization  of  the  medium 
due  to  the  initial  pulse.  From  our  previous  diseussion  we  infer  that  the  actual 
polarization  will  oscillate  at  the  appropriate  plasma  frequency.  In  either  ease 
the  eleetrie  fields  for  t  >  (r^  +  r^  =  R^)/e  are  n°t  at  the  proper  frequency  to  be 
detected  and  we  shall  negleet  the  spurious  terms. 


It  is  eonvenient  to  consider  that  the  incident  beam  intercepts  only  a  finite  seg¬ 
ment  of  the  eylinder.  If  0  denotes  the  angle  between  the  axis  of  the  cylinder  and 
the  direction  of  propagation  of  the  incident  wave,  then  the  length  of  the  cylindrical 
segment  intercepted  will  be  taken  to  be  l  / sin  0  where  i  is  the  width  of  the  beam. 
The  direetly  polarized  wave,  granted  that  it  is  given  by  the  first  Born  approxi¬ 
mation,  will  then  be  stretehed  out  for  a  time  of  the  order  of  (2  H  eot  0)/e.  Thus 

2  l 

i  must  be  ehosen  small  enough  that  the  streteh  time  of  —  eot  0  is  mueh  smaller 
than  2D/(e  sin  0). 


Returning  now  to  the  integral  in  equation  (20),  and  utilizing  all  results  of  the  last 
paragraph  but  one,  then  for  a  fixed  time,  t,  we  have  that  the  integrand  differs 
from  zero  (in  an  essential  way)  in  the  range  of  11^2  as  given  by 


9 


(ri  +  r2)  R12  (rl  +  r2) 

T  +  t  -  -  >  - >  t  -  -  (22) 

c  c  c 

Again  denoting  by  d  the  diameter  of  the  cylinder,  we  find  for 

t  >  2  R  +  (23) 

that  Rl2  >  0  and  the  lower  bound  of  this  important  function  increases  linearly 
with  the  time.  In  the  above  expression,  R  is  the  distance  from  the  receiver  to 
the  closest  point  at  which  the  incident  beam  intercepts  the  cylinder. 


Consider  now  that  wc  hold  the  point  of  first  scattering,  i’p  as  well  as  t  fixed. 
The  points  where  the  second  scatterings  take  place  lie  with  volumes  bounded  by 
the  surface  of  cylinder  and  the  surfaces 


R 


12 


=  T  +  t  - 


(r  !  +  r2)  R12 


t  - 


(i’l  +  r2) 


c  c  c  c 

These  last  two  surfaces  will  be  referred  to  as  the  earlier  and  later  isochronal 

surfaces ,  respectively.  The  later  isochronal  surface  is  actually  defined  only 
2  D 

for  t  >  — — rr  +  2  R  and  it  is  only  for  such  time  that  the  considerations  which  fol- 
sin  0 

low  are  intended. 


(24) 


As  stated  in  the  introduction,  there  is  a  stretching  of  the  pulse  even  if  die  first 
Born  term  alone  is  considered  because  of  the  finite  diameter  of  the  cylinder.  The 
directly  polarized  returns  will  persist  until 


2  R 

t  —  - *  + 


2D 
c  sin  0 


(25) 


2  R  2D 

This  leads  us  to  introduce  die  time  T  -  t - -  — : — rr. 

c  c  sin  0 


In  summary,  for 


T  >  0,  r12  >  0. 
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We  will  be  interested  in  the  dependence  of  the  cross-polarized  field  as  a  func¬ 
tion  of  t/Tq  where  Tq  =  —  ^in  -  .  We  have  that  for  T/Tq  »  1,  R12  ^e  T. 
Also  T/Tq  >  2,  R12  >  c  T0-  Let  us  now  carry  out  the  indicated  differentia¬ 
tions  in  equation  (18).  Note  that  VI  (t  -  Ri2  -  i'i  -  r2)  =  (-  VR12)  DI.  Thus: 


VV 


i_  j_  vr12vr12 


D2  R12  = 


R 


I  -  2  V 


12 


12/ 


V  R  0— I 
12  D 


(26) 


Inserting  these  results  into  equation  (20)  we  get 


ro  [  3a  .a  f  3-* 

£  -  f—  Jd  r2n(r2)Jd  rin<ri> 


A 

(C 
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r12)  (eo  •  R12) 


R 


12 


I  + 


3  -  T 

Rl2D 


R22  D 


-ll 

2 


(27) 


As  T  increases,  die  first  term  in  the  bracket  will  dominate  the  other  terms 
since  the  lower  bound  on  Rj2  increases  linearly  with  T.  This  first  term  gives 
the  radiation  field  resulting  from  the  first  scattering;  the  other  terms  represent 
the  induction  and  quasi-static  fields. 

For  the  time  interval  -c  T+ (r  +r  +  R,  )  <  ct  <  r  +  r  +R,  , 

v  1  2  12;  12  12 


(28  a) 


(28b) 


(28c) 


1  1 


where 


p12 


=  r-i  +  r0  +  R 


12 


in  (28c)  will  be  ignored  as  being 


The  terms  c/ c o  in  (28b)  and  — 
o  v  '  co 
o 

due  to  spurious  d-c  effects  that  are  undetectable  by  the  receiver 


*  -  P12/C 


Inasmuch  as  the  scattered  power  rather  than  the  field  amplitude  itself  is  of  in¬ 
terest  we  shall  eventually  average  over  one  cycle  of  the  r-f  excitation.  Before 
we  do  this  we  shall  perform,  formally,  the  relevant  statistical  averaging. 
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SECTION  3 


STATISTICAL  CONSIDERATIONS 


First  let  us  consider  the  coherent  cross  section.  We  shall  write  n  =  n  +  An, 
where  nQ  is  the  average  particle  density  and  An  is  the  fluctuation.  Statistical 
averages  will  be  denoted  by  the  brackets  \  )> .  Averages  over  a  cycle  will  be  in¬ 
dicated  with  a  bar  over  the  quantity.  The  statistically  averaged  field  amplitude 
is  given  by 


<«>  = 


fr2  r 

-°T-  J 

R 


d3rr 


no  (rl)  no  (r2)  (rl)  A  (r2)^> 


(ec  *  R12^  <go  *  R12^ 
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12 


I  +  — - -  I  + 

R!2  ° 


3  1 

—  I 


R 


2  2 
12  D 


(29) 


If  the  diameter  of  the  cylinder  is  large  compared  with  the  correlation  length,  L„ 
then  <  A(rj)  A  (r2)  >  is  negligible  when  t/t0  »1  for  any  reasonable  corrcla- 

D 

tion  function. 


Squaring  and  averaging  over  one  cycle  of  the  r-f  excitation,  one  has,  for  r/T 
large 
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where  B  (p,  R)  = 


1 

2 


3c 


2  „2 

co  R 


12 


3c 


CO2  R234 


+ 


2 


00 


cos  w(p12-p34) 
(31) 


Each  of  the  terms  in  the  integral  above  can  be  broken  up  into  terms  of  the  form 


J'j  d^r^  d^rg  ^c*  Rjo^  ^12)  x  faction  of  R 42  x  trig  function 


x  |^simi 


of  (COP42)  x  similar  expression  in  R34,  P34 


(32) 


For  t/Tq  sufficiently  large,  R12  (and  R34)  are  bounded  from  below.  Holding  the 
point  of  first  scattering  r^^)  fixed  and  integrating  over  r2  ^4)  will  yield  a  very 
small  result  because  of  the  oscillations  of  the  trigonometric  functions.  This  is 
true  provided  that  n0  does  not  have  a  pronounced  dependence  upon  position.  In 
fact,  we  shall  consider  only  the  case  where  nQ  varies  slowly  over  several  wave¬ 
lengths  of  the  r-f  wave.  However  the  fluctuations  may  vary  considerably  with 
position. 


We  turn  now  to  the  calculation  of  the  incoherent  scattering  and  compute  the 
form  of 


when  t/Tq  >>  1.  In  this  case,  R42  and  R34  >  D  and  if  D/Lc  is  sufficiently 
large,  then  there  is  no  effective  correlation  between  the  density  fluctuations  at 
points  ri  and  r2  and  at  r2  and  Similarly,  with  the  restriction  of  the  initial 
beam  to  a  segment  of  the  cylinder,  and  again  for  t/Tq  large,  correlations 


1 


between  fluctuations  at  17  and  r4  and  at  rg  and  rg  effectively  vanish.  One  has, 
under  these  eircumstanees, 


6 


inc 


74  ^ d3 


r2  d3r3  d3r4  ^ee-:Rl2^ 


(  eo‘  ^34) 


p3  R3 
R  12  R  34 


n  (1)  n  (3)  <A2n>C  (2,4)  +  n  (2)  n  (4)  <A2n>C  (1,3) 


(34) 


+  (<A2n>)2  C  (1,3)  C  (2 , 4) 


B  (R,p) 


2  2 

where  C  (1,3)  is  the  two-point  correlation  function  (we  take  exp  -  R  1Q/L  ). 

Jl  o  c 

In  deriving  this  result  we  took  the  asymptotic  form  of  the  four  point  correla¬ 
tion  function  to  be  a  product  of  two-point  functions.  For  obvious  reasons,  any 
of  the  three-point  correlation  functions  vanish  here. 


Thus,  for  t/T  a  and  D/L  sufficiently  large  the  incoherent  cross  seetion  contains 
only  terms  of  the  Type  I  and  Type  Ille  (in  the  notation  of  Ruffine  and  de  Wolr  b- 
The  terms  of  Type  I,  proportional  to,  say,  n0(l)  nQ(3)  C(2,4),  are  small  for 
reasons  given  in  the  discussion  of  the  coherent  scattering,  namely,  cancella¬ 
tions  due  to  the  trigonometric  functions  in  B  (p,R).  The  same  argument  is 
given  in  (1).  Thus  we  are  led  to  consider  only  terms  of  type  Ille. 


Before  carrying  out  the  integrations  in  equation  (34)  we  turn  to  a  consideration 
of  the  iso-surfaees  mentioned  in  the  previous  seetion. 
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SECTION  4 


THE  ISOCHRONAL  SURFACES 


At  the  position  of  the  scatterer  the  incident  wave  shall  be  considered  as  ap¬ 


proximately  planar  over  the  entire  length  of  the  plasma  cylinder;  we  shall  use 


this  approximation  for  the  rest  of  the  paper.  As  has  been  mentioned  pre¬ 
viously,  the  incident  beam,  for  calculational  purposes,  is  considered  to  be 
spatially  restricted  so  as  to  intercept  the  cylinder  across  a  length  -t/sin  0 
parallel  to  the  cylindrical  axis.  Figure  2  shows  a  cross  section  of  the  cylinder 
in  a  plane  containing  the  axis  of  the  cylinder  and  the  propagation  vector  of  the 
incident  wave.  The  nearest  point  of  interception  of  the  beam  and  cylinder  we 
take  for  the  origin  of  a  cartesian  coordinate  system  £  ,  rj  with  r)  along  the  direc¬ 
tion  of  propagation.  Let  the  point  (£]_,  77^)  represent  a  point  of  first  scattering 
and  (£2*  *7 2 )  ^he  P0*11^  °f  a  possible  second  scattering.  For  the  later  isochronal 
surface,  wc  have 


+  2  R  =  ct 


or 


2  +  <^2  -  ^1> 


2 


(V2  +r)1)  +  C  (T  +  Tq) 


(35) 


(recalling  that  Tfl  =  2D/Csin0). 
0 


Introducing  the  dimensionless  quantities 


;  a  =  —  ;  u 

T0 


V 


(36) 


lfi 


V 


Figure  2.  Geometry  for  Isochronal  Surfaces 


we  get 


1/  = 


l  +  a 
~~2 


2 

u 


2  [a  +  l  -  0] 


0  si 
a  so 


(37) 


The  isosurface  is  a  parabolic  sheet  whose  curvature  is  given  by  -l/(a+  1  -  0) 
and  whose  point  of  interception  with  the  rj  axis  is  ~  cTg,  s*nce’  the 
point  of  first  scattering  0  <  £  <  l ,  the  isochronal  surface  is  actually  a  volume. 
The  distance  t  is  taken  small  enough  that  we  can  ignore  this  complication  in 
what  follows. 


The  upper  and  lower  edges  of  the  cylindrical  cross  section  are  given  by 
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(38) 


vu  =  — +  uctne 


V  =  u  ctn  0 
L 


The  intersections  between  the  parabola  (37)  and  the  lines  (38)  are  at 


uu  =  (a  + 1-0) 


UL  =  (a  +  1-jB) 


ctn  0  dJ  cm  0  + 


a 


a  +  1  -  /3 


2  Q/  +  1 

-  ctn  0  ±4/  ctn  0  +  - - — — 

a  +  1  “  o 


For  a  »  1  -  j3 
± 


u 


U 


(a  +  l  -  3) 


1-8  . 

-  ctn  0  ±  esc  0  -r— 11  sin  0 


u 


±  W  (a  +1-/3)  -  eta  0  ±  CSC  0  ± 


2a 

0 


2  (a  + 1) 


sin  0 


From  this  one  finds 


(39) 


(40) 


UL  -  UU 


sin  0 
.  2 


(41) 


The  distance  between  two  members  of  a  set  of  intercepts,  measured  along  the 
axis  of  the  cylinder  is  approximately  (1  -  cos  0).  If  the  two  intercepts  are 
connected  by  a  chord,  this  chord  makes  an  angle  0  with  the  axis  of  the  cylinder, 
where  0  is  given  by 


1  -  cos  0  ,0^9 

tan  0  =  - : — - —  =  tan  — ;  0  ~ 

sm  0  2  2 


(42) 


since  in  units  of  cTfl  the  diameter  is  l/2  sin  0. 
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The  earlier  isochronal  surfaces  are  easily  obtained  by  replacing  r  by  T  +  T  in 
the  appropriate  equations.  T,  it  mil  be  recalled,  is  the  time  duration  of  the 
initial  r-f  pulse.  Certainly,  when  the  above  approximations  hold  for  later 
isochronal  surfaces  they  hold  for  the  earlier  one.  In  Figures  3,4,  and  5  we 
show  the  later  isochronal  surfaces  for  values  of  /3  =  0,  l/2,  and  1  respectively 
and  for  various  values  of  a. 
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V 


Figure  4.  Isochronal  Surfaces  ( fi  =  1/2) 


V 


CYLINDER 


Figure  5.  Isochronal  Surfaces  (0  =  1) 
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SECTION  5 


THE  INCOHERENT  CROSS  SECTION  FOR  t/T  0  »  1 


In  this  case,  the  dominant  contribution  will  be 


A  A 

where  R^  and  Rg^  are  unit  vectors.  If  the  turbulence  of  the  plasma  is  uniform 
throughout  the  cylinder,  then  the  essential  time  dependence  is  contained  in  the 
product  Rj0  Rg^  occurring  in  the  denominator  of  the  integrand.  This  is  so  be¬ 
cause  the  lower  bound  on  this  product  increases  quadratically  with  the  time; 
thus  the  integral  varies  as  the  inverse  square  of  the  time  for  sufficiently  long 
times.  An  inspection  of  Figures  3-5  also  shows  that  the  scattered,  incoherent, 

cross  polarized  cross  section  becomes  increasingly  aspect  sensitive  as  time  in- 
A  A 

creases  since  Ii^  and  R34  are  asympto  tic  ally  parallel  to  the  cylindrical  axis. 
This  contrasts  with  the  results  in  (1);  in  the  case  treated  in  that  paper  there  is, 
however,  always  an  input. 


These  results  are  what  is  to  be  expected  —  at  the  times  considered  the  radia¬ 
tion  field  from  the  first  scattering  gives  the  largest  contribution  and  the  power 
in  the  radial  field  drops  off  as  the  distance  squared.  Also,  the  polarization  of 
the  first  scattered  wave  is  perpendicular  to  R-^  x  R.  The  amplitude  in  the  R^g 
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direction  is  proportional  to 
the  wave  along  R  onto  e^,  is 


/\  ft 

eQ  x  R12 

/\  ft 
ee  x  R12 


and  the  projection  of  the  polarization  of 
.  When  e0,  ee  and  6j2  lie  in  one  plane 


./\  ft  , 
(eo  x  R12> 


(®e  x  R12> 


<$>  x  ^12)  •  (®e  x  ^12) 


I//\  O  V 
(ee  *  ^12^ 


/A 

(e 


A 

V 


Before  proceeding  we  shall  ehange  our  notation  to  more  nearly  match  that  of 
(1).  First,  at  the  position  of  the  scatterer  the  incident  wave  is  nearly  plane. 

Secondly  a  unit  amplitude  incident  field  is  assumed  in  (1).  We  thus  malic  the 
replacements 


f2 


R 


,  w 
1’  c  P12 
34 


R  '  (r!  +  r2)  +  k  R12 
3  4  34 


— ^  O 

where  k  is  the  direction  of  propagation.  We  multiply  equation  (34)  by  4  n  R  to 
form  a  scattering  eross  seetion.  Finally  it  is  neeessary  to  replaee  4? r  An  by 
the  A  of  (1). 


This  gives  us  a  value  for  the  incoherent  cross  section  when  t/tq  is  not  large 

C7 

and  the  ratio  of  S  =  <  (r  >>  Tn)  >/<  >  for  t/Tq  »  1  and  the  values  of 


<  £  - 
s  me 


>just  diseussed  provides  a  dimensionless  measure  of  the  extent  of 
pulse  stretching. 


_o 

We  have  seen  that  the  incoherent  power  return  should  vary  as  (eT)  when 
t/Tq  >>  1.  Now  if  the  pulse  length,  eT,  is  many  times  larger  than  the 
diameter,  there  will  be  a  time  when  the  scattering  eross  section  should  be 
given  very  nearly  by  the  steady  state  result  of  (1)  Equation  50.  This  will  be 
true  at  times  when  the  later  isosurfaee  has  not  yet  developed  but  the  earlier 
isosurface  is  at  a  large  enough  distance  from  the  initial  point  of  scattering  that 
the  results  of  Ruffine  and  de  Wolf^1)  apply.  (Rceall  that  these  authors  assume 
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that  the  length  of  the  cylinder  is  large  compared  with  the  radius.)  The  "active 
region"  of  the  cylinder  we  are  discussing  is  given  by  the  distance  between  the 
two  earlier  isosurfaces  along  the  axis.  The  only  change  we  must  make  in  equa¬ 
tion  50^)  is  to  replace  the  volume  in  that  formula  by  the  volume  77  a2  -t/sin  Q 
where  a  is  the  radius  of  the  cylinder. 


Making  the  indicated  changes  and  forming  the  ratio  S  we  find 


S  = 


7  t>2_  2 

3(2)  eR  L  sin  6 


(4n)3  77 2  a3  l  L3  ^3  +  sin  $+••••  j 


HP 


d3m 


/\  y\  ys,  /\  S' 

<ec-R12>  <eo'R12>  ^c'  n34*  <e 


eo'R34>J 


R12  R34 


B  (R  1p)  C  (1 , 3)  C  (2 , 4) 


A  A 

b  X  k  o 


Figure  f>.  Geometry  of  Scattering  from  a  Cylinder 


The  relevant  angles  are  shown  in  Figure  6. 

For  the  integration  over  ?3  (r4) ,  we  first  shift  the  origin  to  ri  (r2)  and  use  the 
nature  of  the  correlation  function  to  justify  the  extension  of  the  resulting 
integrals  to  all  of  space.  Furthermore,  we  can  replace  R34  by  R12  except  in 
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the  function  B  (R,  p). 


The  R  „ 
34 


occurring  in  B  (R,  p  )  can  be  approximated  by 


R34  =  R12  +  R12  '  [*4  ‘  ?2  "  3  "  ? 1 >] 

to 

Only  the  terms  in  B  (R,p)  proportional  to  cos  —  (p12  -  o„4)  will  contribute 
(the  one  proportional  to  the  sin  function  gives  zero  in  tills  approximation),  and 
of  these  we  consider  only  terms  independent  of  R-^  and  R34: 


S 


3  (2)7  sin  9,  /_l\  f(  “S  d3f2 

(4TT)3la3  \2/  I'  R2 


(e 


•R 


12' 


<  V  R12) 


2 


1 

^3  +  sin'"'  0  +  •  •  •  | 

The  integral  over  r4  is  essentially  the  r^  volume  (n  a^  L/sin0)  times  the  value 
of  the  integral  over  ?2  when  r-^  is  on  the  axis  of  cylinder.  Thus, 


S 


2  ^  2  /v  /\  /\  a  2 

Assuming  that  a/R12  <<  1,  (ec-R12)  (eo_R12)  ~  (ec_z)  (co‘ z)  or 
cos^  0Q  cos2  0e  in  the  notation  of  (1).  For  the  purposes  of  integration,  we  re¬ 
place  the  isochronal  surfaces  by  planes  perpendicular  to  the  axis  passing 
through  the  points  of  intersection  of  these  surfaces  with  the  axis. 


Then 


S 


6a 


[cos2  0q  cos^  0c] 


dz 

~2 

z 


Left-hand 
Iso- volume 

(Equation 


Right-hand  (44) 

Iso-volumc 

continued  on  next  page) 
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(44)  (cont. ) 


3a  cos  0  cos  0c 
 o  


Left 

hand 


+ 


hand 


2  2 

3a  cos  0Q  cos  0c 


2  c  T 


(c  T  +  D 
l  sin  0  11 


cT+cT +  .  _ 

sin  0 


Thus  the  pulse  stretching  falls  off  as  the  inverse  square  of  the  stretch  time. 

VI.  NON-UNIFORM  ELECTRON  DENSITY 

In  order  to  obtain  Equation  (44)  a  number  of  simplifying  assumptions  were  made 
including: 


o  2 

1)  That  the  correlation  function  is  of  the  form  exp  (-R  / 1  c) 

2 

2)  That  the  fluctuation  strength  <  A  >  is  a  constant. 


Since  neither  of  these  conditions  are  observed  in  actual  range  observations,  it 
is  interesting  to  see  what  results  when  they  arc  dropped.  We  shall  also  ask  a 
somewhat  different  question,  namely,  given  that  at  some  time,  t,  wc  observe 
a  direct  and  a  cross  polarized  return,  what  is  the  error  in  assuming  that  the 
two  scatterings  which  caused  the  cross-polarized  return  both  occurred  in  the 
volume  that  caused  the  direct  return?  We  assume  that  the  incident  pulse  is 
much  longer  than  the  diameter  of  the  cylinder  and  is  infinitely  wide.  The 
direct  scattering  takes  place  in  the  illuminated  volume  of  the  cylinder  defined 
as  the  volume  lying  between  planes  located  at  distances  c  (t  +  T)/2  and  ct/2. 
When  the  radar  line  of  sight  is  along  the  cylinder  areas,  at  least  one  of  the 
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collisions  must  occur  in  this  volume.  Since  we  need  only  be  concerned  with 
the  case  of  small  angle-0,  we  take  this  to  be  true  in  what  follows.  Proceeding 
as  before,  Equation  (43)  can  be  put  into  the  form 


-  r 

2  _  l_  _o 

’  inc  '  2  2 

R 


(45) 


A  a  2  a  2 

F(k-1^)  x  (e^k)  (e^k)  , 


k  =  kR12,  F  (k)  C  (R)  cos  (k  •  R)  d  R. 


The  limits  on  r^  and  r2  can  be  derived  from  (21).  When  r^  and  r2  lie  in  the 
illuminated  volume,  (45)  is  equivalent  to  the  result  obtained  in  (1).  (The 
factor  of  l/2  arises  from  the  fact  that  we  have  always  taken  r^  and  ro  to  be  in 
the  same  correlation  cell. )  We  now  proceed  to  estimate  the  contribution  to  (45) 
when  r^  lies  outside  of  the  correlation  volume. 

For  sufficiently  large  separation  between  scatterers,  k  lies  along  the  cylinder 
axis.  Defining  this  axis  by  a  unit  vector,  b,  we  have 


F  (kc  +k  b)  F  (1^  -  kb)  (ec-b)2  (e^b)2 


(4G) 


/ 


dIl  % 


V'2  R12 


<A(1)2^  <A(2 ))> 


(45) 
(cont. ) 


where  is  the  volume  of  the  cylinder  below  the  plane  at  ct/2  and  V2  is  the  il¬ 
luminated  volume.  To  estimate  the  value  of  the  second  integral  in  (46),  choose 
the  origin  of  the  coordinate  system  at  the  center  of  the  illuminated  volume  and 


replace  R  by  r^.  The  integral  becomes 


// 


d£i  di2 


R‘ 
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12 


v2  <A(2)2^>  tt  a2 
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■lr 


c  T 


4  cos  0 


<A(12)^ 


dxi 


where  iQ  is  the  distance  from  the  start  of  the  cylinder  to  the  center  of  the  il¬ 
luminated  volume.  Thus,  we  may  define  a  new  quantity"  S'  to  be  the  ratio  be¬ 
tween  the  return  in  which  one  scattering  is  outside  the  illuminated  volume  to 
the  return  in  which  both  scatterings  are  outside  the  illuminated  volume. 


S'  = 


37 T2 

a 

f  <  A“  >  dx 

16 

<A(2)2>  j 

1  x2 

/\  A  A  A2  A  AO 

F  Gs0  +kb)  F  (j^j  -  kb)  (eo-b)  (e^b)2 
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d  k  F  (kQ  +k)  F  (1^-k) 


(y  jj  frg 

sin  0 


where  cos  0  =  k  •  b  . 
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For  the  special  case  of  a  Gaussian  correlation  function,  S'  reduces  to 

a 


S'  =  6 


2  2 
cos  0  cos  0 
o_ c 


2  o  2  2 

3  +  cos  0„  +  cos  0^+9  cos  0O  cos  0C 


o 


<  A  (2)*  > 


c  T 


4  cos  0 


<  A  (x)  >  dx 
2 


(Please  see  Reference  1  for  more  detail. ) 


Finally,  we  assume  that  both  unit  polarization  vectors  have  the  same  projection 
on  the  plane  containing  1^  and  b  and  we  estimate  the  contribution  from  the  adjoin¬ 
ing  resolution  cell.  Assuming 


cos  1  jio  '  b  «  20 
cos  0Q  =  cos  @2  =  .22 
c  T 

f o  2  cos  0 


then 


S' 


2  x  (.  0025) 
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c  T 

2  cos  0 


<A(i)2> 

<A  (2)2  > 


Taking  a/ 1  «*0. 1  (surely  an  overestimate) 


«s  5  x  10 


<A(12)  > 
<A(2)2  > 


we  have  S' 


9  2 

The  ratio  of  <  A  (1)^  >/<  A(2)  >  can  be  estimated  from  the  ratio  of  direct 
polarized  wake  scattering  in  adjacent  resolution  cells.  Thus,  if  S'  <  10-2  for 
example,  the  cross-polarized  scattering  can  be  assumed  to  arise  from  the 
same  resolution  cell  as  the  directly-polarized  scattering.  This  will  occur 
when 


A(1)2>  <5X102 
A(2)2> 

or  when  the  direct  polarized  scattering  decreases  by  less  than  approximately 
23  db  per  resolution  cell. 
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